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Abstract 


The  design  of  heat  exchangers  traditionally  focuses  on  the  known  con¬ 
straints  of  the  problem  such  as  inlet  and  outlet  temperatures,  flow  rates,  and 
pressure  drops.  This  leads  mainly  to  a  sizing  problem  where  the  designer 
must  select  surfaces,  flow  configuration,  and  materials  to  meet  the  mini¬ 
mum  design  objectives.  An  alternate  approach  based  on  an  acceptable  level 
of  thermodynamic  irreversibility  (entropy  generation)  has  been  proposed. 
When  the  entropy  generation  level  hats  been  set,  the  geometric  parameters 
of  the  heat  exchanger  can  be  determined.  The  design  of  a  plate-fin  type, 
gas-to-gais  recuperator  for  a  regenerative  open  Brayton  cycle  has  been  used 
as  a  demonstrative  device.  The  resulting  heat  exchanger  designs  are  then 
examined  to  determine  what  caused  the  differences  and  why  either  method 
should  be  preferred  over  the  other. 
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Chapter  1 


INTRODUCTION 


1.1  Background 

Energy  conservation  is  a  topic  that  has  recently  received  considerable  atten¬ 
tion.  Available  work  has  been  recognized  as  a  valuable  commodity  and  its 
destruction  in  seemingly  otherwise  efficient  engineering  processes  has  also 
seen  a  renewed  increase  in  interest.  The  use  of  second  law  analysis  and  ther¬ 
modynamic  irreversibility  minimization  have  been  proposed  as  techniques 
that  should  become  an  integral  part  of  the  design  of  engineering  processes 
and  components  [ij.  It  is  not  apparent  what  the  relationship  of  traditional 
methods  of  design  and  these  irreversibility  minimization  techniques  will  be. 

1.2  Conventional  Design  Procedures 

The  design  of  a  component,  such  as  a  heat  exchanger,  generally  involves  the 
sizing  of  that  component  to  meet  specified  performance  parameters  within 
known  constraints.  The  heat  exchanger  area  and  volume  are  usually  de¬ 
signed  to  be  the  minimum  required  to  meet  the  specifications  as  this  will 
also  usually  be  the  best  design  from  an  ecomonic  aspect.  The  fluid  inlet  and 
outlet  temperatures,  flow  rates,  and  pressure  drops  are  usually  specified  and 
it  is  the  designer’s  task  to  determine  construction  type,  flow  arrangement, 
materials  and  surfaces  to  meet  those  requirements  [2],  The  effect  of  the 
component  design  on  the  overall  system  performance  should  also  be  consid- 


ered  to  ensure  the  tradeoffs  performed  within  its  design  are  still  valid  when 
integrated  into  the  system. 


1.3  Irreversibility  Minimization 

Any  heat  transfer  process  is  generally  accompanied  by  thermodynamic  ir¬ 
reversibility  or  entropy  generation.  The  entropy  generated  is  in  direct  pro¬ 
portion  to  the  amount  of  useful  work  dissipated  in  that  process.  If  this 
irreversibility  can  be  minimized  within  a  particular  component  of  a  power 
cycle,  the  useful  power  output  of  the  cycle  should  increase  [3l.  Bejan  has 
shown  how  the  entropy  generation  rate  can  be  reduced  in  a  counterfiow  gas- 
to-gas  heat  exchanger  Al,  and  it  is  this  method  that  will  be  used  to  create 
designs  that  will  be  compared  to  conventionally  designed  heat  exchangers, 
and  their  impact  on  overall  cycle  performance. 
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Chapter  2 


SELECTION  OF 
BRAYTON  CYCLE 
PARAMETERS 

2.1  Cycle  Components 


A  regenerative  open  Brayton  cycle  was  chosen  as  the  vehicle  to  test  the 
different  design  methods  for  the  recuperator.  A  schematic  of  the  cycle  com¬ 
ponents  is  shown  in  Figure  2.1. 


Regenerative  Cycle  Schematic 

Figure  2.1 


A  compressor  polytropic  efficiency  was  assumed  and  used  throughout 
all  calculations,  T}cp  =  0.95.  Similarly,  a  turbine  polytropic  efficiency  was 
assumed,  =  0.90.  For  ease  of  calculations,  only  one  fluid  was  used, 
i.e.,  air  was  modeled  as  a  perfect  gas  with  Cp  =  1.0  kj/kg  K  and  R  = 
0.287  kj/kg— °  K  both  assumed  constant.  The  combustor,  then,  was  treated 
as  a  perfect  heat  transfer  device,  but  was  not  an  injection  point  for  fuel. 

2.2  Cycle  Parameters 

A  reasonable  compressor  inlet  temperature  was  selected,  Ti  =  300°  if,  along 
with  an  inlet  pressure  of  Pi  =  1.126  x  10^  iV/m^.  A  nominal  value  of  turbine 
inlet  temperature  was  also  selected,  T4  =  1300° if.  The  selection  of  an 
appropriate  pressure  ratio  required  the  following  consideration. 

The  addition  of  a  recuperator  to  a  simple  open  Brayton  cycle  causes  a 
shift  in  the  pressure  ratio  for  maximum  cycle  efficiency  from  that  of  the  non- 
regenerative  cycle.  A  heat  exchanger  effectiveness  of  £  =  0.8  was  selected 
and  cycle  efficiency  without  losses  was  calculated  for  pressure  ratios  varying 
from  4  to  10.  A  sample  calculation  is  contained  in  Appendix  A,  results  are 
plotted  in  Figure  2.2. 


C'/C:e  Ef'^iciency  vs  Pressure  Rati< 

Hgur#  2.2 


S'  0.485  ^  / 
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Based  on  these  results,  a  pressure  ratio  of  P-i/ P\  =  5  was  selected.  With 
this  pressure  ratio  determined,  the  maximum  cycle  efficiency  without  losses, 
rjoo)  could  be  calculated.  This  quantity  will  be  used  for  comparison  purposes 
and  was  computed  as  rjoo  =  0.5760. 

2.3  Heat  Exchanger  Characteristics 

A  plate-fin  heat  exchanger  surface  was  selected  from  those  cataloged  in  ref¬ 
erence  [5].  Its  characteristics  and  dimensions  are  shown  in  Figure  2.3.  It  has 
an  uncomplicated  geometry  and  a  reasonably  wide  range  of  Reynolds  num¬ 
bers  in  which  consistent  results  can  be  expected.  The  material  is  aluminum, 
with  a  thermal  conductivity  of  kf  =  190  w/m  K. 


F»n  pitch  *  2  0  p€r  lO  ■  78.74  p«r  m 
Plate  spacing,  b  ■  0.750  m  »  19.06  *  l0"^rT> 

Fin  length  •  12.0  in  -  304.8  x  lO'^m 

Flow  passage  hydraulic  diameter,  4r^  =  0.0474  ft  »  14.453  x  10'^m 

Fin  metal  thickness  *  0.032  in.  aluminum  =  0.813  x  10‘^m 

Total  heat  transfer  area/volume  between  plates,  3  *  76.1  ft*/ft^  *  249  672  m^/m^ 

Fin  ar99JxoX9\  area  ’  0.608 


Plain  Plate-Fin  Sur-^ace  2.0 

Figure  8.3 
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CONVENTIONAL 
DESIGN  METHOD 


3.1  Introduction 

This  chapter  outlines  the  conventional  design  procedure  that  was  used  to 
size  the  single-pass,  counterflow  gas-to-gas  heat  exchanger.  This  procedure 
was  performed  for  various  combinations  of  heat  exchanger  effectiveness  and 
specified  pressure  drops.  A  complete  numerical  example  is  provided  in  Ap¬ 
pendix  B.  Final  results  for  this  conventional  method  are  shown  in  Table 
3.1. 


3.2  Conventional  Methodology 

For  a  single-pass  counterflow  design,  when  the  core  dimensions  on  one  side 
are  fixed,  the  dimensions  on  the  other  side  are  also  known.  This  means  the 
design  is  driven  by  the  side  that  has  the  more  stringent  AP/P  requirement. 
The  method  of  determining  the  controlling  side  is  contained  in  Appendix 
B.  In  every  cjwe  considered,  the  hot  side  of  the  heat  exchanger  was  the 
controlling  side  and  so  its  requirements  drove  the  design. 

The  following  is  a  step-by-step  heat  exchanger  sizing  design  procedure 
which  closely  follows  that  outlined  in  reference  [2]. 


1.  For  the  given  heat  exchanger  effectiveness,  determine  the  fluid  outlet 


temperatures.  Calculate  the  fluid  mean  temperature  on  each  side  and 
evaluate  fluid  physical  properties  />,,  po,  [\l p)m,  and  p. 


2.  Determine  Ntu  for  the  exchanger  and  then  Nm  for  each  side.  The  in¬ 
fluence  of  longitudinal  heat  conduction  is  ignored  in  this  first  iteration 
of  design. 

3.  Estimate  hot  and  cold  side  pressure  drops,  select  an  appropriate  Nrc, 
and  then  a  value  of  j//  from  Figure  2.3. 

4.  Calculate  mass  velocity,  G,  from  information  in  steps  1-3  and  the  cor¬ 
responding  value  of  AP/P. 

5.  Calculate  Nr^  and  determine  values  of  j  and  f  from  Figure  2.3. 

6.  Compute  heat  transfer  coefficient,  h;  temperature  effectiveness  of  the 
fins,  fjy;  and  the  total  surface  temperature  effectiveness,  r\a. 

7.  Calculate  heat  transfer  tirea.  A;  minimum  free  flow  area,  Ao',  heat  ex¬ 
changer  frontal  area,  A/r;  flow  length,  L;  and  heat  exchanger  volume, 
Vol. 

8.  Compute  AP/P  from  known  conditions  and  calculated  parameters. 

3.3  Analysis 

For  each  case  show  in  Table  3.1,  cycle  efficiency  with  losses,  r],  was  calculated 
and  the  cycle  efficiency  degrade,  listed  for  each  case.  A  correlation 

between  heat  exchanger  effectiveness,  e,  total  allowed  pressure  drop,  k,  cycle 
efficiency  degrade,  and  heat  exchanger  volume  was  sought. 

The  relationship  between  cycle  efficiency  degrade  and  effectiveness  at 
various  levels  of  k  is  shown  in  Figure  3.1.  The  increased  level  of  degrade 
at  lower  values  of  effectiveness  is  as  expected,  and  for  a  particular  value  of 
effectiveness,  the  degrade  increases  with  increasing  pressure  drops. 

The  effects  of  varying  effectiveness  and  pressure  drops  on  heat  exchanger 
volume  is  depicted  in  Figure  3.2.  The  higher  the  effectiveness,  the  longer 
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and  therefore,  larger  the  heat  exchanger.  Similarly,  the  greater  pressure 
drop  for  a  particular  effectiveness  also  increases  the  volume. 

Figure  3.3  was  derived  by  extracting  values  of  effectiveness  at  nominal 
values  of  rjoo  —  V  plotting  them  on  the  appropriate  k  curve  of  Figure  3.2. 
This  resulted  in  curves  of  constant  ~  V  that  have  a  minimum  volume 
at  a  particular  value  of  k  and  e.  It  is  emphasized  that  these  “minimum” 
volumes  were  obtained  at  constant  values  of  cycle  degrade  with  varying 
pressure  drops  and  effectiveness.  These  “minimum”  volumes  are  plotted 
against  effectiveness  in  Figure  3.4  with  the  respective  k  values  annotated  on 
the  graph. 


Case 

6 

k 

Voo  V 

L(m) 

Vol  (m^) 

C-1 

1-0 

0.02 

0.0043 

— 

— 

C-2 

1.0 

0.04 

0.0089 

— 

— 

C-3 

1.0 

0.06 

0.0137 

— 

— 

C-4 

1.0 

0.08 

0.0188 

— 

— 

C-5 

0.9 

0.0 

0.0449 

— 

— 

C-6 

0.9 

0.0240 

0.0498 

13.69 

3.69 

C-7 

0.9 

0.0388 

0.0550 

13.82 

2.64 

C-8 

0.9 

0.0558 

0.0603 

14.04 

2.19 

C-9 

0.9 

0.0739 

0.0660 

14.65 

1.98 

C-10 

0.8 

0.0 

0.0833 

— 

— 

C-11 

0.8 

0.0192 

0.0886 

6.15 

1.12 

C-12 

0.8 

0.0367 

0.0941 

6.56 

0.85 

C-13 

0.8 

0.0550 

0.0999 

6.72 

0.71 

C-14 

0.8 

0.0740 

0.1060 

7.15 

0.65 

C-15 

0.7 

0.0 

0.1165 

— 

— 

C-16 

0.7 

0.0188 

0.1221 

3.75 

0.53 

C-17 

0.7 

0.0365 

0.1278 

3.95 

0.39 

C-18 

0.7 

0.0565 

0.1338 

4.17 

0.34 

C-19 

0.7 

0.0734 

0.1401 

4.32 

0.31 

C-20 

0.6 

0.0 

0.1455 

— 

— 

C-21 

0.6 

0.0186 

0.1512 

2.52 

0.29 

C-22 

0.6 

0.0376 

0.1571 

2.69 

0.22 

C-23 

0,6 

0.0557 

0.1632 

2.85 

0.19 

C-24 

0.6 

0.0764 

0.1695 

2.98 

0.17 

Table  3.1:  Conventional  Method  Results 
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Chapter  4 

MINIMUM  ENTROPY 
DESIGN  METHOD 

4.1  Introduction 

The  use  of  irreversibility  concepts  in  the  design  of  the  gas-to-gas  counterflow 
heat  exchanger  will  be  presented  in  this  chapter.  For  the  same  mass  veloci¬ 
ties  that  were  used  in  Chapter  3,  the  heat  exchanger  geometric  parameters 
can  be  determined  to  arrive  at  minimum  levels  of  entropy  generation.  The 
outputs  of  this  design  procedure  are  heat  exchanger  effectiveness  and  the 
pressure  loss  suffered  on  each  side.  For  a  more  detailed  explanation  of  this 
procedure,  the  reader  is  directed  to  references  [4]  and  [6j. 

4.2  Background 

Since  the  method  deals  with  the  entropy  generation  within  the  recuperator, 
it  is  appropriate  to  begin  with  an  entropy  flux  analysis  on  the  control  volume 
of  the  recuperator  in  Figure  4.1.  The  heat  transfer  from  the  outer  walls  is 
assumed  to  be  negligible.  The  entropy  generation  rate  can  then  be  written 
as 

•S'  —  5fn)l  ~  ^in)2 


' 


‘^trean  1 


A  A 


Stream  2 


Recuperator  Control  Volune 

Figure  4.1 


Substituting  expressions  for  the  ideal  gas  entropy  changes  based  on  pressure 
and  temperature  relationships,  this  becomes 


(4.2) 


5  =  Cm,nMTiout/Tl)  +  iR/Cp)iln{Pi/Piout)] 

+Crr^zMT2out/Tj)  +  {R/C,)2ln(P2/ Piout)] 

The  number  of  entropy  production  units,  N,,  is  defined  as 

N.  =  S/Cmaz  (4.3) 

Writing  expressions  for  the  first  law  of  thermodynamics  of  the  control  vol- 

(4.4) 


ume. 


Cmm(7’l  -  Tiout)  +  C'm,„(T2  -  Tjout)  =  0 

and  heat  exchanger  effectiveness, 

e  =  [Tiout  -  Ti)/{T2  -  Ti) 
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(4.5) 


permits  the  elimination  oiTiout  and  T2oiit  in  equation  (4.2)  which  allows  the 
rate  of  entropy  production  formula  to  be  nondimensionalized  as 


Kv? 


N,  - 


(4.6) 


A  special  form  of  this  last  expression  can  be  obtained  in  the  case  of  nearly 
ideal  heat  exchangers.  When  the  stream-to-stream  AT’s  and  frictional  A  P’s 
are  held  small,  two  inequalities  can  be  applied. 


«  1 

(4.7) 

<<  1 

(4.8) 

Applying  these  and  expressing  c  in  terms  of  Ntu  yields 


^  ^min  ,  Pj  ,  , 
N,  cs  — - in—  +  In 


1  - 


laz  V  Pz  /  . 


1  - 


(1  -  T 

exp  -  iV| 

(1  - 

\  max  / 

1  -  ( 

^mas  \ 

1  —  ^f^exp 

-A(u  (1  -  ' 

.  \  m  a  1  /  J 

*  fc  ■  « 


*■; 


R 


AP' 

~P 


R 

a 


AP^ 


(4.9) 


1  '  -  '  *  \'"P/  2 

The  first  two  terms  can  be  characterized  as  a  contribution  from  capacity 
rate  imbalance,  the  third  represents  contribution  due  to  finite  Ntu,  and  the 
last  two  show  fluid  friction  effects. 

Requiring  the  streams  to  be  balanced  and  applying  the  calculus  of  limits 
as  Cmin  Cmaz,  equation  (4.9)  can  be  reduced  to 


C 

A-'  rn.n.T. 


In'^  +  In 
Ti 


1-^ 

Cmaz  V  T2 


iTi-  iTi 
V  V  T'. 


Cmaz  \ 


AP\ 


When  the  overall  Nt^  is  expressed  in  terms  of  Ntu  for  each  side  of  the  heat 
transfer  surface, 

J_^J_ 

iVtu  Ntui  Cjnax  Ntu2  '  ' 

the  number  of  entropy  generation  units  can  be  divided  into  three  contribu- 


N,  —  N ,  imbalance  "b  ^sl  "b  ^a2 


(4.12) 


Applying  the  limit  Cmin  ^max>  the  imbalance  terms  of  equation  (4.10) 
can  be  written  as 


imbalance  — 


0(1- 


,  ^min  ,  T'l 
Cmax  Pi 


(4.13) 


which  will  clearly  vanish  when  Cmin  —  Cmax-  Since  this  is  the  case  for  the 
design  here,  the  number  of  entropy  production  units  for  each  side  can  be 
determined  from 


N,2=  \ 


Tx  V  ^2 ;  Nt^x 

®'j_ 

Pi  V  P2  /  ^tu2 


AP\ 

P  )x 


(4.14) 


AP\ 


4.3  Design  Considerations 

The  number  of  entropy  generation  units  that  have  been  deemed  acceptable 
for  a  particular  design  is  a  function  of  the  hydraulic  radius,  minimum  free 
flow  area  and  flow  length.  It  is  necessary,  then,  to  express  N,  in  terms  of 
these  and  other  known  flow  parameters.  Using  the  definitions  of  Ntu  and 
friction  factor  for  each  side 

^tuX,2  =  (b/r;,)iV,(  (4.16) 


(AP/P)i.2  -  f{L/r^]G^/{2pP) 
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in  equations  (4.14)  and  (4.15)  leads  to 


where 


+f>Bf{L/n)g^ 

[L/rh)N,t 


ai  =1 


(4.18) 


1^2  —  ^min/ 


(4.19) 


^>1  — 


62  =  1 


(4.20) 


B,  =  (i?/Cp)i 


Bt  =  (B/Cp): 


VV  Tx  \Ti) 


(4.21) 


(4.22) 


g  ^  G/(2pP)J 


(4.23) 


As  can  be  seen  in  equation  (4.18),  the  ratio  L/r^^  performs  a  trade-off  func¬ 
tion,  i.e.,  for  a  fixed  g  and  Nftf,  there  will  be  an  optimum  L/r^,  which  results 
in  a  minimum  Nsi,2-  When  the  function  N,i  2  is  minimized,  the 
is  given  by 

[L/rKUt  =  G  (^/(l/p.P,  +  \/phPh^ 

4.4  Minimum  Entropy  Design  Methodology 

The  method  used  s  based  on  optimizing  (L/r^).  Since  r^^  is  set  by  the  selec¬ 
tion  of  the  plate-fi:'  surface,  this  yields  a  length  that  will  produce  minimum 
entropy  for  the  given  hydraulic  radius. 


The  following  outlines  the  methodology  employed.  A  complete  numerical 
example  is  shown  in  Appendix  C.  Final  results  for  all  cases  are  tabulated  in 
Table  4.1. 

1.  For  the  cases  listed  in  Table  3.1  that  resulted  in  a  heat  exchanger 
volum.e,  use  calculated  G  and  to  determine  {L|T^^)opt■ 

2.  Compute  resultant  heat  exchanger  effectiveness  and  hot  and  cold  side 
pressure  drops. 

3.  Calculate  N,^,  and 

4.  Calculate  the  heat  exchanger  volume  for  each  case. 

5.  Determine  the  cycle  efficiency  degrade,  fjoo  ~ 

4.5  Analysis 

Figure  4.2  shows  the  nearly  linear  relationship  of  cycle  efficiency  degrade  and 
heat  exchanger  effectiveness.  There  is  not  the  clear  definition  of  different 
values  of  pressure  drop  ais  was  observed  for  the  conventional  method  in 
Figure  3.1.  In  fact,  there  is  even  a  wider  range  of  k  values  here,  from 
k  =  0.056  for  case  E-6  to  k  =  0.269  for  case  E-24,  and  yet  the  results  are 
almost  linear  across  that  range. 

The  volume  plotted  in  Figure  4.3  is  that  volume  obtained  by  minimizing 
entropy  generation  in  the  heat  exchanger  and  should  not  be  interpreted  as 
a  minimum  volume,  but  rather  as  an  “optimum”  volume  for  minimizing 
entropy  generation. 

To  achieve  a  better  basis  for  comparison  of  the  two  methods,  the  entropy 
generated  by  the  various  designs  must  be  discussed. 


€ 


Case  #  Vol(m^)  k 


^oo 


E-6 

32.03 

8.63 

0.0563 

0.955 

0.0347 

0.0325 

E-7 

25.86 

4.93 

0.0727 

0.945 

0.0440 

0.0417 

E-8 

22.04 

3.44 

0.0876 

0.935 

0.0528 

0.0503 

E-9 

20.12 

2.72 

0.1014 

0.926 

0.0606 

0.0583 

E-11 

24.35 

4.46 

0.07658 

0.942 

0.0461 

0.0441 

E-12 

19.08 

2.46 

0.1069 

0.922 

0.0641 

0.0614 

E-13 

16.11 

1.07 

0.1320 

0.908 

0.0784 

0.0760 

E-14 

15.08 

1.38 

0.1559 

0.896 

0.0917 

0.0897 

E-16 

19.80 

2.77 

0.996 

0.927 

0.0596 

0.0573 

E-17 

15.40 

1.53 

0.1422 

0.902 

0.0844 

0.0823 

E-18 

13.19 

1.07 

0.1778 

0.884 

0.105 

0.1028 

E-19 

12.11 

0.86 

0.2058 

0.872 

0.1211 

0.1189 

E-21 

16.98 

1.93 

0.1252 

0.913 

0.0741 

0.0718 

E-22 

12.87 

1.04 

0.1799 

0.882 

0.1067 

0.1048 

E-23 

11.52 

0.76 

0.2251 

0.862 

0.1335 

0.1309 

E-24 

10.47 

0.60 

0.2685 

0.847 

0.1596 

0.1541 

Table  4.1:  Entropy  Generation  Method  Results 


Chapter  5 

COMPARISONS  AND 
RESULTS 

5.1  Introduction 

To  this  point,  the  minimum  entropy  method  results  have  been  compared  to 
the  conventional  method  results  in  terms  of  conventional  method  standards, 
i.e.,  pressure  drops,  efficiencies,  and  effectiveness.  It  is  necessary  to  compare 
these  methods  in  terms  of  the  minimum  entropy  method,  i.e.,  the  entropy 
generated  in  each  design. 

5.2  Conventional  Method  Entropy  Generation 

The  number  of  entropy  generation  units  for  each  design  can  be  computed 
from  previously  determined  information  and  equations  (4.14)  and  (4.15). 
An  example  of  these  calculations  is  contained  in  Appendix  D.  Results  for 
the  conventional  cases  are  listed  in  Table  5.1.  The  large  deviation  in  N,  and 
loo  —  1  for  cases  C-21  through  C-24  are  a  result  of  violating  the  assumption 
that  1  —  e  <<  1.  These  conventional  cases  will  be  discarded  from  further 
consideration.  However,  the  minimum  entropy  cases  at  corresponding  mass 
velocities  are  still  valid. 

When  the  entropy  generation  units  are  plotted  against  the  degrade  in 
cycle  efficiency,  Figure  5.1,  the  dependence  on  pressure  drop  is  still  present. 


The  grouping  of  data  points  is  caused  by  the  discrete  changes  in  heat  ex¬ 
changer  effectiveness,  decreasing  from  left  to  right. 

Figure  5.2  was  created  by  using  the  same  values  of  rjoo  —  V  that  were  used 
for  Figure  3.4.  The  value  of  N,  at  each  of  these  cycle  efficiency  degrades 
was  read  off  at  a  particular  value  of  k.  Plotting  these  “minimum”  volumes 
against  entropy  generation  rate  exhibits  the  general  trend  that  the  larger 
the  volume  of  the  heat  exchanger,  the  lower  the  total  entropy  generated. 
This  is  in  consonance  with  the  results  predicted  by  reference  4^ 

5.3  Minimum  Entropy  Generation 

The  cycle  efficiency  degrade  vs.  entropy  generated  for  the  minimum  entropy 
method  is  similar  in  overall  shape  to  that  of  the  conventional  method,  but, 
as  Figure  5.3  shows,  the  relationship  is  nearly  linear  and  does  not  show  the 
discrete  pressure  drop  differences. 

Figure  5.4  shows  that  the  “optimum”  volumes  determined  hy{L/r^^)oJ,^ 
follow  essentially  the  same  path  as  that  of  the  conventional  method.  For  a 
direct  comparison,  Figures  5.2  and  5.4  are  plotted  together  on  Figure  5.5. 
To  ensure  that  the  {L/rh)opt  expression,  (4.24),  was  giving  minimum 


.V/ 


■»> 


entropy  results,  a  series  of  calculations  was  performed  at  {L/r^^)  values  on 
both  sides  of  the  optimum  for  case  E-9.  Figure  5.6  shows  that  the  minimum 
entropy  solution  was,  in  fact,  being  determined.  It  is  also  noted  that  the 
magnitudes  of  N,  and  r/oo  —  t]  are  essentially  the  same  over  this  range  of 
heat  exchanger  effectiveness. 


5.4  Comparisons 


The  calculations  performed  on  case  E-9  to  prove  minimization  also  yielded 
heat  exchanger  volumes  at  a  constant  mass  velocity  G.  These  additional 
data  points  are  plotted  on  Figure  5.7.  The  dash-dot  line  connecting  them 
is  a  locus  of  recuperator  volumes  along  a  constant  G  path.  This  shows  why 
the  volumes  at  minimum  IV,  or  -  V  along  a  constant  G  path  are  larger 
that  the  miminum  volume  of  a  given  N,  or  rj^  -  rj. 
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.A   ^  1. 


Figure  5.8  shows  a  comparison  between  conventional  "minimum’’  vol¬ 
umes  and  minimum  entropy  “optimum”  volumes.  The  minimum  entropy 
method  yields  volumes  that  are  larger  than  the  conventional  method  for  the 
same  cycle  efficiency  degrade. 
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Case  # 


C-2 

0.0089 

— 

C-3 

0.0137 

— 

C-4 

0.0188 

— 

C-5 

0.0449 

— 

C-6 

0.0498 

0.0445 

C-7 

0.0550 

0.0494 

C-8 

0.0603 

0.0550 

C-9 

0.0660 

0.0610 

C-10 

0.0833 

— 

C-11 

0.0886 

0.0900 

C-12 

0.0941 

0.0967 

C-13 

0.0999 

0.1036 

C-14 

0.1060 

0.1108 

C-15 

0.1165 

— 

C-16 

0.1221 

0.1505 

C-17 

0.1278 

0.1584 

C-18 

0.1338 

0.1669 

C-19 

0.1401 

0.1748 

C-20 

0.1455 

— 

C-21 

0.1512 

0.2308 

C-22 

0.1571 

0.2405 

C-23 

0.1632 

0.2501 

C-24 

0.1695 

0.2607 

Table  5.1:  ( 

IJonventional  Method  Entropy  Generation 
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Figure  5.5 
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Chapter  6 


CONCLUSIONS 


6.1  Discussion 

The  nature  of  the  results  and  conclusions  drawn  for  the  example  cycle  would 
be  the  same  for  other  cycles  and  recuperator  surfaces.  For  high  effective¬ 
ness  recuperators,  the  counterflow  arrangement  is  dictated.  Here  with  the 
same  finned  surface  on  both  sides,  the  (AP/P)c/(AP/P)/,  ~  (P/i/P^)*  or 
approximately  1/25.  If  a  closer  plate  spacing  were  used  on  the  clean  air  cold 
side  more  AP/P  would  appear  on  the  cold  side  for  any  total  EP/P.  This 
would  result  in  even  smaller  minimum  volumes  at  any  given  N,  or  rice  ~  V- 

6.2  Conclusions 

In  selecting  a  gas  turbine  recuperator,  essentially  the  same  results  are  ob¬ 
tained  either  by  the  conventional  method  of  minimizing  degradation  of  cycle 
efficiency,  rjeo  —  r],  or  by  minimizing  entropy  generation,  N,. 

The  selection  of  an  “optimum”  recuperator  requires  that  care  be  exer¬ 
cised  in  the  selection  of  the  path  along  which  the  optimization  is  calculated. 
It  has  been  shown  that  the  volume  of  the  recuperator  for  minimum  N,  or 
Hoo  -  1  along  a  constant  mass  velocity  path  is  much  greater  than  the  true 
minimum  volume  at  any  ri^o  —  rj  oi  N,. 
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Appendix  A 

Sample  Calculation  for 
Effect  of  Recuperator  on 
Cycle  Efficiency 


Ts  =  — -  =  820.1°fc 

1300 


Compressor: 


Pj  \  inpc) 

A/ 


Recuperator: 


Cycle  Efficiency: 


/  P2  \  iilc)h) 


=  514.2°fc 


T3-T2 

'  Ts  -  Tj 

T3  =  (€)(T5  -  T:)  +  Tj 


Tz  =  (.8)(820.1  -  514.2) 514.2  =  695, 2°A: 


(74  -  Ts)  +  (Ti  -  T2)  (1300  -  820.1)  -r  (300  -  514.2) 


74-T3 


1300  -  695.2 


=  0.4651 


WWW 


u 


Appendix  B 

Conventional  Design 
Method  Numerical  Example 


Before  calculations  can  be  performed,  two  points  need  to  be  clarified:  1)  the 
method  of  determining  how  loss  of  turbine  work  due  to  pressure  drops  in 
the  system  is  distributed  and;  2)  the  method  for  determining  which  side  of 
the  heat  exchanger  has  the  more  stringent  AP/P  requirement. 

Turbine  work  is  given  by  the  expression 


tht  =  mrjtCpTH  1  “ 


mrjtCpTh 


Pc  +  APc 


Ph-APk 


PcV\l+^ 


Expanding  the  last  term. 


APi,  APc 


so  that  loss  of  turbine  work  can  be  expressed  as 

PcV  f^Ph  APc 


‘•iV, 


r.wwr.' 


and  it  is  clear  from  Figure  B.l  that  determines  the  loss  of  turbine 

work  irregardless  of  how  it  is  distributed.  Introducing  the  definitions  r  =  ^ 


specific  entropy 

Temperature  -  Entropy  Diagram 
F igure  B.l 


and  9  —  (^)  ,  it  can  be  shown  that  for  a  counterflow  exchanger  where 

\  ■‘«  /  avg 


Ac  =  Ah  and  Lc  =  Lh, 


APh/Ph  _  9 


APc/Pc  r2 

Using  a  slightly  different  form  of  the  last  turbine  work  expression, 


Athf 

mrjtCpTh 


'APh/Ph  1  A  Pc 
-APc/Pc  \  Pc 


Awt  ^  A  ^ 

mrjtCpTh  Vr^  /  Pc 


Introducing  =  k,  expressions  for  hot  and  cold  side  pressure 

drops  are, 

_  k[e/r^) 

Ph  1  +  ^/r2 

AP,  _  k 

P,  ~  TTJp 

The  numbered  steps  that  follow  correspond  to  those  outlined  in  Section 
3.2.  The  input  design  parameters  used  for  the  example  calculations  were 
e  =  0.8  and  k  =  0.04. 

B.l  Temperatures  and  Fluid  Properties 


Ti  =  300°  ifc 

T2  _  f 

“  \KJ 

so  that 


T2  =  300(5)  <  =  486.8°it 


Ts  = 


(^) 


npt 


(^) 


T4 


Ts 


(4.8)(-^)(^) 


T  -1 


1300 


=  868.5°fc 


so  that 


Similarly 


e  = 


T3-T2 


n-Ti 

T3  =  e{Ts-T2)  +  T2 

T3  =  (.8)(868.5  -  486.8)  =  792.2°^: 

Te  =  Ts-eCTs-Tz) 


Te  =  868.5  -  (.8){868.5  -  486.8)  =  563. Tfc 

Since  this  is  a  counterflow  design  with  C*  >  0.5,  the  best  choice  for 
average  temperatures  is  the  arithmetic  average  on  each  side 


Tem  — 


T2  +  Tz  486.8  +  792.2 


=  639.5°A: 


2  2 

Inlet,  outlet  and  average  densities  are  next  computed.  Pv  =  RT  or 
V  =  RT / P  but  V  =  1/p,  so  p  =  P/ RT 

1.126  X  lO^V/m^ 


PS 


{287 j/kg  ifc)(868.5°A:) 


=  0.4517fc^/m'’ 


50 


1.126  X  10^ 


(287)(563.1) 


=  0.6967 kg /m^ 


P'  m.A 


=  i  (—+  —'j  =  1.824m VA:^ 
2  \  P5  P6  / 


5.63  X  10® 


(287)(486.8) 


=  4.030A:g/ m® 


5.63  X  10® 


(287)(792.2) 


=  2.476fcg/m® 


=  0.3260m®/it5 


B,2  Nt 


For  a  couterflow  heat  exchanger, 


j  _  1  -  exp[-jVt^(l  -  c‘)] 

1  -  c*eip(-A^tu(l  -  c')] 


but  for  the  balanced  flow  case  (c*  =  1)  under  design, 


1  +  ^(u 


.r  _  f  _  0-8  _  . 

^fu  —  -  —  -  —  4 

1  -  c  1-0.8 


Since  gas  is  the  working  fluid  on  both  sides,  it  was  estimated  that  both  sides 
would  have  approximately  the  same  surface  resistance  so  that 


Ntu,h  =  ^«u.c  -  2^■tu  =  8 


B.3  Estimated  Pressure  Drops 


Recall 


1.1193 


avg 
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0.03862 


APh  _  k{e/r^)  _  (0.04)(27.983) 
~P^  ~  1  +  "  1  +  27.983 


APc  _  k  0.04 

~Pr  ~  TVffr^  ~  1  +  27.983 


0.00138 


From  Figure  2.3,  an  estimate  of  j//  =  0.0035/0.0077  was  obtained 


B.4  Mass  Velocities 

The  mass  velocity  on  each  side  is  estimated  from  the  relationship 

G=\(  (^£IL)  " 

As  a  first  approximation,  rjo  was  assumed  to  be  0.90. 

G  -  [/'(2)(1126x  10^iV/m^)(0.9)\  /0.03862\  /0.0035 

^  A  (1.824m3/fc^)(0.7)J  j  V  8  /  VO. 0077 

G^^  =  n  .586kg  /  —  s 


^  _  rA2)(5.63  X  10^)(0.9)  W0.00138W0.0035\1  = 

'  “  [v  (.3260)(0.7)l  J  V  8  /  V 0.0077/ 

Gc  —  n  .582kg/m^  —  s 

G^Dh  (17. 583A:j/m^  —  s)(0.0474/t)(. 3048m  ft) 

3.396  X  \Q~^kg/m  -  s 


,V„  =  (1:  583)(.0474)(.3048)  ^ 

3.159x10-5 

Entering  Figure  2.3  with  these  values  yields 

fh  =  0.0078 


jh  =  0.00355 


f,  =  0.0077 


jc  =  0.0035 


B.5  Heat  Transfer  Coefficients  and  Fin  Effective¬ 
nesses 

With  the  Colburn  factor  known,  the  heat  transfer  coefficient  can  be  deter¬ 
mined  from 

h  =  jGcpP7'^ 

This  result  is  used  to  compute  the  fin  characteristic  length, 


The  temperature  effectiveness  of  the  fin  is  next  calculated  from 

tanhml 


A  good  approximation  for  I  is  given  by  half  the  plate  spacing  minus  the  fin 
thickness. 

The  total  surface  temperature  effectiveness  can  then  be  calculated 

r?o  =  1  -  (1  -  r]f)Aj/A 

where  A;j A  is  a  characteristic  of  the  surface  selected. 

(0.00355)(17.586fc(7/m2  -  s){\<^kijkg  k) 

hh  = - 2 - 

(0.7)  J 
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ing  resistances  and  wall  thermal  resistances, 


(0.9848)(79.19)  (.9850)(78.06) 


=  0.02583 


=  38.714u;/m-  k 


B.6  Dimensions 


The  heat  transfer  surface  area  is  computed  Irom 


A  = 


Nt.C 


A  =  (^)(36Q0%Mr)(l-Q^>/fcg  fe)  ^  ,03 

38.714M;/m^  k 

Minimum  free  flow  area  is  given  by 

A,  =  (W/G)  =  „  =  0.05686m> 

°  ^  ^  I7.586kg/m^  -  s 

In  order  to  determine  frontal  area,  <t  must  be  computed  from  the  relationship 

d  -l-  d  -f*  2t 

Here  a  parting  plate  thickness  of  0.4mm  was  assumed,  so 

(0.75m)(76.1/t7/t3)(0.0474/4/0 
0.75»n  +  0.75m  +  2(.4mm)(. 1/2.54  in/mm) 


a  =  0.44162 


A  It  —  Ac/ 17 


0.05686m*  , 

A,t  =  -  =  0.12876m^ 

^  0.44162 


Appendix  C 

Minimum  Entropy  Design 
Method  Numerical  Example 


This  example  uses  the  data  from  case  C-12  of  Table  3.1, 

G  =  \2A^&kg/rn^  -  s,  j  =  0.00355,  and  /  =  0.0078. 

C.l  Optimum  Length 

In  order  to  use  equation  (4.24),  intermediate  values  must  first  be  determined 


j  0.00355 
"  (iVp,)2/3  =  (5^  "  ° 

The  temperature  span  parameter,  r,  cannot  be  calculated  until  Ts  is  known, 
so  an  estimate  based  on  previous  iterations  is  used.  If  the  calculated  value 
of  r  is  close  to  the  value  used,  the  computations  are  valid. 


G  \  Bf{i/p,P,  +  l/p^P^) 


0.365/0.00450 


(0.287)(0.0078)  (5-1^ - 


0.325  ,  1  818 

63x  10^  ^  1  126x  in' 


Since  Th,  is  set  by  the  geometry  selected,  L  can  be  calculated 

=  5283(.0474/<)(.3048m/7i:) 


Lapt  =  19.08m 

C.2  Heat  Exchanger  Effectiveness  and  Pressure 
Drops 


From  equation  (4.16), 


-  {L/r^)N,t 


Ntu.ti  =  (5283)  (0.00450)  =  23.79 
Applying  the  same  eissumption  as  was  used  in  Chapter  3, 


Ntu  =  23.79/2  =  11.89 

Then,  since  this  is  still  a  balanced  counterflow  design, 

Ntu  11.89 

1  + Ntu  ~  1  +  11.89  “ 

For  hot  and  cold  side  pressure  drops,  equation  (4.17)  is  used. 
AP 

(-^)^  =  f{Wrk)9V{2pP)  ^  f{LlrH)9^ 


AP 


=  (0.0078)(5283)(0.05005)^  =  0.010322 


=  (0.0077)(5283)(0.00946)^  =  0.00369 


=  0.10322  +  0.00369  =  0.1069 


C.3  Entropy  Generation 


:  '-n 


\ 


BftTiVif. 


Equations  (4.14)  and  (4.15)  will  be  used  to  determine  entropy  generation  on 
both  sides. 


Nsh 


Nsh  --- 


^tu,h 

0.365 


R 


+  - 


AP' 
~P  yK 


23.79 


+  (0.287)(0. 10322) 


Nsh  =  0.04497 


0.365 


23.79 


+  (0.287)(0.00369) 


Ns.  =  0.0164 


Ns  =  Nsh  +  Nse 


Ns  =  0.04497  +  0.0164  =  0.0614 


C.4  Heat  Exchanger  Volume 


Since  Afr  is  calculated  from  minimum  free  Row  area  and  other  surface  geo¬ 
metrical  constraints, Ayr  from  the  corresponding  case  in  Chapter  3  was  used 
to  calculate  the  heat  exchanger  volume. 


Vol  =  L{Afr)  =  19.08m(0.12876m^)  =  2.46m^ 


C.5  Cycle  Efficiency  Degrade 


Temperatures  that  reflect  the  actual  effectiveness  and  pressure  drop  must 
first  be  determined. 


Ps 


Pz 


y{\  -  k)  =  5(1  -  0.1069)  =  4.465 


60 


(4.465)  1 

1300 


=  884.8°fc 


Ts  —  €(Ts  —  Ti)  -r  T2 


Tz  =  0.922(884.8  -  486.8)  +  486.8 


Tz  =  853.8°jt 


(T4  -  Ts)  +  (Ti  -  Tz) 


T4-T3 


1300  -  884.8)  +  (300  -  486.8)  „ 

n  =  i i =  0.5119 
'  1300  -  853.8 


T]^-r]  =  0.5760-  0.5119  =  0.06412 


A  check  on  the  validity  of  the  t  value  assumed  shows 


Tz  [Tz 

Tz  V  Tz 


486.8  V  884.8 


r  =  0.367 


Therefore,  assumed  value  of  0.365  was  valid 


Appendix  D 

Sample  Calculation  of 
Entropy  Generation  for 
Conventional  Designs 


Example  uses  case  C-12.  From  previous  example, 

r  =  0.34461 


APfc 

Pk 


0.0359 


AP. 

Pc 


0.00136 


^tu,h  —  ^tu.e  —  8 

Using  equations  (4.14)  and  (4.15), 


N,,h  =  - - +  (0.287)(0.0359)  =  0.05338 

8 

Similarly, 

+  (0.287)(0.00136)  =  0,04347 

8 


N,  =  =  0.09685 

Adding  the  first  term  of  each  side’s  expression  gives  the  AT  contribution, 
and  the  sum  of  the  second  terms  is  the  AP  contribution  to  the  entropy 
generation. 


Ns^AT 

=  0.08615 

N,,ap 

=  0.01070 
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*0 


h> 
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